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We  present  quality  factors  of  single-defect  photonic-crystal  resonant  cavities  with  asymmetric  cladding  layers. 
The  resonators  studied  here  are  dielectric  slabs  patterned  with  two-dimensional  photonic  crystals  on  a  sap¬ 
phire  substrate.  Three-dimensional  finite-element  and  finite-difference  time-domain  routines  were  used  to 
analyze  the  electromagnetic  properties  of  these  cavities.  We  observe  that  high  quality  factors  (—800)  can  be 
obtained  in  these  cavities  for  reasonable  structures  with  thick  enough  dielectric  slabs.  This  work  was  moti¬ 
vated  by  the  need  to  place  photonic-crystal  resonators  on  a  substrate  to  improve  heat  dissipation  in  photonic- 
crystal  lasers.  ©  2002  Optical  Society  of  America 
OCIS  codes:  140.0140. 


Photonic-crystal  microcavity  lasers  are  potentially  attrac¬ 
tive  optical  sources  in  future  communication  systems. 
They  have  several  advantages  such  as  lithographically 
defined  wavelengths1  and  expected  low  operating  powers. 
Much  work  remains  to  be  done,  however,  for  these  sources 
to  find  mainstream  applications.  Single-defect  photonic- 
crystal  lasers  were  first  demonstrated  in  pulsed  mode  at 
low  temperatures.2  Room-temperature  pulsed  operation 
has  since  been  demonstrated  in  lasers  with  larger  mode 
volumes,1,3-6  and  continuous  wave  (cw)  room- 
temperature  operation  has  been  reported  in  a  photonic- 
crystal  laser  with  a  fusion-bonded  substrate  that  had  a 
layer  that  was  subsequently  converted  to  Al,.Oy  ,T  The 
cw  demonstration  used  a  large  hexagonal  cavity  formed 
in  a  triangular  lattice.  Continuous  wave  operation  of 
single-defect  photonic-crystal  lasers,  however,  has  not 
been  reported  to  date.  To  operate  such  a  laser  at  room 
temperature,  it  is  necessary  to  design  a  high-Q  (quality 
factor)  resonant  cavity  that  dissipates  heat  well.  Poor 
heat  dissipation  is  likely  the  major  factor  preventing  cw 
operation  of  single-defect  lasers  in  undercut  device  geom¬ 
etries.  As  in  the  cw  operation  of  larger  photonic-crystal 
cavities,  a  strategy  for  improving  heat  dissipation  in 
single  defect  cavities  is  to  wafer  bond  the  semiconductor 
membrane  containing  the  two-dimensional  photonic  crys¬ 
tal  to  a  low-index-of-refraction,  high-thermal-conductivity 
substrate.  Here  we  choose  sapphire  as  the  substrate  ma¬ 
terial.  This  complicates  the  electromagnetic  design, 
however.  There  are  two  major  effects  of  adding  a  uni¬ 
form  substrate  below  the  photonic-crystal  membrane. 
First,  the  presence  of  a  substrate  breaks  the  symmetry  of 
the  structure  about  the  midplane  of  the  device.  This  re¬ 
moves  the  ability  to  classify  modes  as  either  even  or  odd 
and  removes  the  bandgap  in  the  guided  modes  of  the 
membrane.8  Second,  the  substrate  also  will  increase  op¬ 
tical  losses  in  the  cavity  because  of  the  larger  radiation 
loss  into  the  substrate  compared  with  an  air-clad  struc¬ 
ture.  Here  we  investigate  the  quality  factor  of  these 
resonators  as  a  function  of  the  semiconductor  membrane 


thickness.  We  will  show  that  it  is  possible  to  obtain 
high-Q  resonant  cavities  by  increasing  the  membrane 
thickness  while  maintaining  the  in-plane  mode  localiza¬ 
tion.  A  significant  complication  occurs,  however,  when 
the  slab  becomes  thick  enough  to  support  multiple  modes. 
The  multiple  slab  modes  close  the  bandgap  in  the  guided 
modes.8,9  This  introduces  radiation  loss  in  the  in-plane 
direction.  We  analyze  single-layer  and  multiple-layer  di¬ 
electric  slabs  as  a  function  of  the  slab  thickness  and  com¬ 
pare  the  calculated  threshold  material  gain  required  to 
achieve  lasing  in  a  single-defect  cavity  for  each  case.  As 
an  aside,  we  note  that  it  should  also  be  possible  to  obtain 
high-Q  resonant  cavities  in  the  presence  of  a  substrate  if 
the  in-plane  mode  localization  is  relaxed,  but  this  issue  is 
not  pursued  here. 

Electromagnetic  modes  in  purely  two-dimensional  pho¬ 
tonic  crystals  can  be  classified  as  having  either  TE  or  TM 
polarization.  TE  modes  have  the  electric  field  polarized 
in  the  plane  of  periodicity  and  TM  modes  have  the  mag¬ 
netic  field  polarized  in  the  plane  of  periodicity.  Once  the 
photonic  crystal  becomes  finite  in  height,  this  classifica¬ 
tion  scheme  is  lost.  For  example,  in  the  original 
photonic-crystal  laser  demonstration,2  the  modes  were 
not  purely  TE  or  TM.  That  demonstration,  however,  had 
inversion  symmetry  about  the  midplane  of  the  device  be¬ 
cause  the  cladding  layers  were  symmetric.  This  symme¬ 
try  allows  a  classification  of  the  modes  as  being  even  or 
odd  about  the  midplane.  At  the  midplane,  the  even  mode 
is  TE  and  the  odd  mode  is  TM,  but  away  from  that  plane 
this  TE  or  TM  classification  is  no  longer  valid.  In  a  cav¬ 
ity  with  symmetric  cladding  layers  there  is  a  gap  formed 
for  the  even  guided  modes  as  well  as  for  the  odd  guided 
modes.  The  bandgaps  in  the  even  and  odd  mode  spec¬ 
trum  do  not  completely  overlap,  but  since  the  two  sets  of 
modes  are  orthogonal  there  is  still  a  gap  for  each  mode 
classification.  A  sapphire  substrate  breaks  the  symme¬ 
try  and  there  is  no  way  to  classify  the  modes  of  the  cavity. 
As  a  result  there  is  no  longer  a  gap  in  the  guided  mode 
spectrum  of  the  dielectric  membrane.  This  has  the  po- 
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tential  to  limit  or  even  eliminate  the  ability  to  obtain  in¬ 
plane  localization  for  a  mode  formed  by  a  defect  in  the 
two-dimensional  photonic-crystal  lattice. 

A  more  severe  problem  associated  with  the  sapphire 
cladding  is  increased  radiation  loss  into  the  sapphire. 
The  radiation  loss  out  of  the  plane  of  the  photonic  crystal 
is  predicted  to  be  the  dominant  optical  loss  in  the  demon¬ 
strated  single-defect  photonic-crystal-membrane  laser 
structures.  This  loss  arises  because  part  of  the  optical 
mode  is  not  totally  internally  reflected  at  the 
semiconductor-air  interfaces  at  the  top  and  bottom  of  the 
resonant  cavity.  In  other  words,  the  localization  of  the 
mode  in  real  space  causes  a  spread  of  the  mode  in  wave- 
vector  space,  and  a  portion  of  the  mode  lies  above  the 
light  line  in  k  space.9-12  This  means  that  the  mode  is  not 
completely  guided  by  the  membrane.  This  problem  will 
become  worse  in  a  structure  bonded  to  a  sapphire  sub¬ 
strate  because  the  light  line  of  the  substrate  has  a  lesser 
slope  than  the  light  line  in  an  air-clad  structure.  This  re¬ 
duces  the  separation  between  the  mode  in  k  space,  which 
is  centered  about  the  M  point  and  the  light  line. 

Figure  1  shows  the  geometry  of  the  resonant  cavity. 
Figure  2  shows  the  calculated  band  structure  for  both  the 
symmetric  air-clad  and  the  asymmetric  air-sapphire-clad 
photonic-crystal  membranes.  These  calculations  were 
performed  by  a  finite-element  analysis.  Whitney  1  forms 
were  used  as  vector  basis  functions  in  the  finite-element 
calculation.  The  mesh  density  used  was  2948  tetrahe¬ 
dral  elements  per  cubic  lattice  constant.  The  symmetric 
air-clad  structure  consisted  of  a  three-layer  dielectric 
membrane  surrounded  by  air.  The  membrane  core  had 
an  index  of  refraction  of  3.4  and  a  thickness  of  0.2a, 
where  a  is  the  lattice  constant.  The  dielectric  cladding 
layers  on  both  sides  of  the  core  had  a  thickness  of  0.1a 
and  an  index  of  3.1.  The  ratio  r/a  was  0.3,  where  r  is  the 
radius  of  the  holes.  The  air-sapphire-clad  structure  had 
the  same  dielectric  membrane  structure  but  with  a  sap¬ 
phire  layer  below  the  membrane  with  an  index  of  refrac¬ 
tion  of  1.6.  The  dielectric  slab  supports  only  a  single 
mode  at  this  thickness.  The  calculations  illustrate  the 
effects  on  the  radiation  mode  condition  when  a  substrate 
is  used  as  a  cladding  layer.  A  larger  optical  loss  results 
for  the  same  membrane  and  photonic-crystal  lattice  com¬ 
pared  with  the  symmetric  air-clad  geometry.  There  is 
also  a  shift  in  the  bands  from  the  air-clad  structure  to  the 
sapphire  substrate  structure  and  a  reduction  in  the  band- 


Fig.  1.  Illustration  of  resonant  cavity  geometry  modeled  in  this 
work. 


Fig.  2.  Finite-element  calculations  of  the  dispersion  relations  of 
the  photonic-crystal  membrane  for  (a)  Symmetric  air  cladding, 
(b)  asymmetric  air-sapphire  cladding. 

gap.  This  occurs  because  for  the  same  membrane  thick¬ 
ness,  the  modes  in  the  sapphire-bonded  structure  have  a 
larger  effective  index.  This  results  in  a  slight  lowering  of 
the  mode  frequency.  The  air  band  has  a  larger  frequency 
reduction  than  the  dielectric  band,  resulting  in  the  gap  re¬ 
duction. 

Figure  3  shows  a  cross-sectional  view  of  the  electric 
field  intensity  inside  a  photonic-crystal  cavity  bonded  to 
sapphire.  The  figure  was  obtained  from  a  three- 
dimensional  finite-difference  time-domain  (FDTD)  simu¬ 
lation.  The  FDTD  calculations  used  a  rectangular  mesh 
with  25  points  per  lattice  constant.  The  figure  clearly 
shows  an  increase  in  the  radiation  loss  into  the  substrate 
compared  to  the  radiation  loss  into  the  air  above.  This 


Kim  et  al. 


Vol.  19,  No.  8/August  2002/J.  Opt.  Soc.  Am.  B  1779 


increased  loss  will  significantly  reduce  the  Q  of  the  reso¬ 
nant  cavity. 

The  quality  factors  of  the  resonator  were  calculated  by 
using  the  FDTD  code.10'13-15  Here  we  are  considering  a 
triangular  lattice  with  a  single  missing  hole,  and  it  is  the 
Q  of  the  doubly-degenerate  donor-resonant  mode  associ¬ 
ated  with  this  missing  hole  that  we  are  concerned  with. 
The  Q  of  the  cavity  was  obtained  by  multiplying  the  reso¬ 
nant  frequency,  which  is  obtained  from  a  Fourier  trans¬ 
form  of  the  field,  by  the  energy  stored  in  the  cavity  and 
then  dividing  by  the  power  lost  out  of  the  cavity.  Follow¬ 
ing  Ref.  10,  we  can  separate  this  Q  into  a  Q  due  to  radia¬ 
tion  loss  out  the  top  of  the  cavity,  Q  lop ;  radiation  loss  into 
the  sapphire,  Q  bottom  ;  and  radiation  loss  in  plane  out  of 
the  photonic  crystal,  Qs ide.  In  our  calculations,  the 
power  radiated  in  plane  was  calculated  by  integrating  the 
normal  component  of  the  Poynting  vector  over  the  four 
sidewalls  of  the  calculation  domain.  The  height  of  these 
walls  was  5  d/a,  where  d  is  the  dielectric  slab  thickness 
and  a  is  again  the  photonic  lattice  constant.  The  angle  at 
which  power  radiated  was  divided  into  vertical  radiation, 
and  in-plane  radiation  therefore  varies  slightly  with  dla . 
This  does  not  significantly  affect  the  results,  however. 
The  reciprocal  total  Q  is  just  the  sum  of  the  reciprocal 
component  Qs.  This  total  Q  was  verified  by  comparing 
the  result  of  the  Q  obtained  from  the  time  rate  of  change 
of  the  energy  in  the  resonant  cavity. 

Figure  4  shows  the  calculated  Qsidc ,  Q  bottom  >  and  total 
Q  as  a  function  of  the  slab  thickness.  As  noted  above, 
r/a  =  0.3,  and  seven  lattice  periods  on  each  side  of  the  de¬ 
fect  were  included.  Figure  4(a)  shows  the  data  for  a 
three-layer  dielectric  slab  on  sapphire,  and  Fig.  4(b) 
shows  the  data  for  a  single-layer  dielectric  slab  with  a  re¬ 
fractive  index  of  3.4  on  sapphire.  The  membrane  that 
was  considered  in  Fig.  4(a)  was  an  InGaAsP  core  of  refrac¬ 
tive  index  3.4  with  InP  cladding  on  both  sides,  which  had 
an  index  of  refraction  of  3.2.  The  figure  shows  that  large 
Q  values  can  be  obtained  for  thick  membranes.  As  ex¬ 
pected,  Q bottom  increases  as  the  slab  thickness  increases. 
Q  bottom  will  be  limited  in  practice  by  the  maximum  obtain¬ 
able  etch  depth.  It  should  also  be  expected  that  Qbottom 
depends  on  r/a.  Figure  5  shows  the  calculated  Qtop  and 
Qbottom  f°r  a  fixed  membrane  thickness  of  d/a  =  1  as  a 
function  of  r/a .  The  corresponding  values  for  a  symmet¬ 
ric  air-clad  structure  are  included  for  comparison.  The  Q 
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Fig.  3.  FDTD  calculation  of  electric  field  intensity  inside  the 
air-sapphire-clad  photonic-crystal  resonant  cavity. 


(a) 


(b) 

Fig.  4.  FDTD  calculation  of  Qbottom ,  Qside ,  and  Qtotal  as  a  func¬ 
tion  of  the  membrane  thickness  for  (a)  three-layer  slab  on  sap¬ 
phire,  (b)  single-layer  slab  on  sapphire. 


r/a 

Fig.  5.  FDTD  calculation  of  Qtop ,  Qverticai  and  Qbottom  as  a  func¬ 
tion  of  r/a  with  the  membrane  thickness  fixed  at  1.6  a. 

values  for  the  symmetric  structure  included  in  Fig.  5  are 
for  radiation  loss  into  the  top  (or  bottom,  but  not  both)  in 
order  to  simplify  the  comparison  with  the  sapphire- 
bonded  structure.  As  expected,  the  radiation  loss  into 
the  sapphire  is  larger  than  the  radiation  loss  into  the  air 
for  both  the  symmetric  air-clad  and  sapphire-bonded  reso¬ 
nant  cavities.  Figures  4(a)  and  4(b)  also  show  that  for 
thick  dielectric  slabs  the  total  Q  of  the  cavity  is  not  domi¬ 
nated  by  Qbottom  at  all,  but  by  Qside  •  T°  obtain  the  maxi¬ 
mum  value  of  Q  in  the  cavity,  there  is  therefore  an  impor¬ 
tant  trade-off  between  Qbottom  an<i  Qside-  This  occurs 
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because  as  the  thickness  of  the  dielectric  slab  increases,  it 
supports  more  modes.  These  modes  close  what  would  be 
a  gap  in  the  guided  modes  of  the  slab  if  the  slab  were  sym¬ 
metric.  From  Fig.  4(a)  and  4(b)  we  see  that  Qside  decays 
slightly  more  rapidly  with  increasing  slab  thickness  for 
the  three-layer  slab  than  it  does  for  the  single-layer  slab. 
This  is  true  because  the  higher-order  slab  modes  are 
spaced  farther  apart  in  frequency  for  the  three-layer  slab 
than  they  are  for  the  single-layer  slab.  This  has  the  ef¬ 
fect  of  closing  the  gap  more  quickly. 

Figures  4  and  5  illustrate  the  fact  that  good  vertical 
confinement  can  be  obtained  in  sapphire-bonded 
photonic-crystal  resonant  cavities  by  increasing  the  mem¬ 
brane  thickness  and  choosing  the  r/a  ratio  carefully.  We 
also  found  that  Qside  is  a  strong  function  of  the  slab  thick¬ 
ness,  so  d/a  must  be  chosen  carefully. 

We  also  checked  that  the  in-plane  confinement  was  not 
affected  by  the  asymmetric  cladding  of  the  sapphire.  We 
plot  in  Fig.  6  the  Qside,  Qtop,  Qbottom,  and  total  Q  values 
as  a  function  of  the  number  of  photonic-crystal  lattice  pe¬ 
riods  for  a  membrane  thickness  of  0.85a  and  an  r/a  value 
of  0.3.  This  membrane  supports  a  single  mode.  The  fig¬ 
ure  shows  that  the  Qside  values  increase  rapidly  as  the 
number  of  lattice  periods  increases,  until  the  total  Q  is 
dominated  by  the  out-of-plane  losses.  It  is  important  to 
remember,  however,  that  Qside  becomes  limited  owing  to 
the  closing  of  the  bandgap  as  d/a  increases.  This  closing 
of  the  bandgap  with  increasing  d/a  is  not  a  result  of  the 
asymmetry  introduced  by  the  sapphire  layer,  and  we  have 
therefore  chosen  to  consider  conditions  in  Fig.  6  that  il¬ 
lustrate  the  effects  of  the  asymmetry  more  simply.  Since 
this  in-plane  loss  does  not  saturate  in  Fig.  6  with  increas¬ 
ing  number  of  lattice  periods,  there  appears  to  be  no  evi¬ 
dence  of  significant  coupling  between  the  confined  reso¬ 
nant  mode  and  the  band  edge  states  that  correspond  to 
those  bands  that  have  odd  symmetry  in  the  air-clad  ge¬ 
ometry.  We  expect  that  as  the  membrane  thickness  in¬ 
creases,  the  influence  of  the  sapphire  cladding  on  the 
resonant  electromagnetic  mode  decreases,  since  the  thick¬ 
ness  of  the  InP  cladding  layer  between  the  InGaAsP 
waveguide  core  and  the  sapphire  layer  is  increasing.  As 
the  InP  layer  thickness  increases  to  infinity,  the  modes 
become  even  and  odd  modes  of  the  InGaAsP- InP  semi¬ 
conductor  layer  structure. 

Finally,  we  calculated  the  threshold  material  gain  re¬ 
quired  for  these  single  defect  resonant  cavities.  As 


Number  of  Photonic  -  Crystal  Layers 

Fig.  6.  FDTD  calculation  of  Q  values  as  a  function  of  the  num¬ 
ber  of  photonic-crystal  periods  for  a  single-mode  slab. 


Fig.  7.  FDTD  calculation  of  threshold  modal  gain  as  a  function 
of  d/a  for  the  three-layer  and  single-layer  slab  cases. 


shown  in  Figs.  4(a)  and  4(b),  the  single-layer  slab  has  a 
larger  total  Q  value  than  the  three-layer  slab  for  the  same 
slab  thickness.  A  laser  design  must  also  account  for  the 
confinement  factor,  however,  which  is  larger  for  the  three- 
layer  dielectric  slab  than  it  is  for  the  single-layer  struc¬ 
ture.  Figure  7  shows  the  calculated  threshold  material 
gain  versus  d/a  for  each  of  the  two  slab  configurations. 
This  was  obtained  using  the  relation 


Q  ~  “fph  “ 


2  7m  eff  1 

k  r^th 


where  T  is  the  confinement  factor,  neff  is  the  effective 
modal  index,  and  gth  is  the  threshold  material  gain.  This 
model  considers  only  radiation  losses,  of  course.  T  was 
calculated  for  the  resonant  mode  by  using  the  FDTD  code. 
The  vertical  confinement  factor  was  calculated  over  a  re¬ 
gion  of  40  nm  at  the  center  of  the  dielectric  slab.  40  nm 
corresponds  to  two  grid  points  at  this  mesh  density.  The 
in-plane  confinement  factor  was  taken  to  be  approxi¬ 
mately  1.  We  see  that  there  is  a  tradeoff  between  Q  and 
the  confinement  factor  and  that  the  lowest  threshold  gain 
condition  occurs  for  a  single-layer  slab  at  d/a  =  1.4. 
This  is  not  necessarily  the  maximum  Q  value,  but  the 
maximum  value  of  the  product  T  Q .  We  note  that  a  ma¬ 
terial  gain  value  of  2200  cm-1  for  four  quantum  well’s,  is 
a  reasonable  value  for  this  material  system.16  At  d/a 
—  1.4,  finite-difference  calculations  indicate  that  the  di¬ 
electric  slab  supports  two  modes:  the  fundamental  mode 
and  a  second  mode  that  has  approximate  odd  symmetry 
about  the  midpoint  of  the  slab.  The  radiative  coupling 
into  this  higher-order  mode  will  be  very  small  in  a  laser 
with  a  quantum  well  active  region  near  the  center  of  the 
slab  because  the  confinement  factor  for  this  mode  is  very 
small. 

We  show  that  high-Q  single-defect  photonic-crystal 
resonant  cavities  can  be  formed  by  bonding  a  semiconduc¬ 
tor  membrane  to  a  high-thermal-conductivity  sapphire 
substrate.  The  increased  radiation  loss  into  the  sapphire 
can  be  compensated  for  by  increasing  the  semiconductor 
membrane  thickness.  There  are  trade-offs,  however,  be¬ 
tween  Qbottom  and  Qs ide  as  the  slab  thickness  increases. 
We  see  no  effects  of  the  asymmetry  in  the  in-plane  Q.  Fi¬ 
nally,  there  is  an  additional  trade-off  between  the  confine¬ 
ment  factor  and  the  total  Q  value  in  designing  a  laser  cav¬ 
ity. 
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